Sum rule for response function in nonequilibrium Langevin systems 
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We derive general properties of the linear response functions of nonequilibrium steady states 
in Langevin systems. These correspond to extension of the results which were recently found in 
Hamiltonian systems [A. Shimizu and T. Yuge, J. Phys. Soc. Jpn. 79, 013002 (2010)]. We discuss 
one of the properties, the sum rule for the response function, in particular detail. We show that 
the sum rule for the response function of the velocity holds in the underdamped case, whereas it 
is violated in the overdamped case. This implies that the overdamped Langevin models should be 
used with great care. We also investigate the relation of the sum rule to an equality on the energy 
dissipation in nonequilibrium Langevin systems, which was derived by Harada and Sasa. 

PACS numbers: 05.70.Ln, 05.40.-a 



I. INTRODUCTION 

Responses to weak perturbations contain important in- 
formation on the physical properties of systems, and are 
often used to characterize their states. General prop- 
erties of the response function are useful in investigat- 
ing the consistency of obtained results both in theories 
and in experiments, and in making prediction from them. 
For equilibrium states many general properties are well- 
known, which are based on the linear response theory [l[ . 
Much interest has been devoted recently to extending the 
properties to nonequilibrium steady states (NESSs) |2j- 
fl2j . Several properties of the linear response function of 
NESSs were found by Shimizu and the author in gen- 
eral Hamiltonian systems [l2j]. These properties consist 
of experimentally measurable quantities only, which is 
a distinguished feature from the formal results in many 
other works. One of the general properties is the sum 
rule for the response function (l2j : 



Re$i(«;F)^ = ({B,A}p) 



F,0' 



(1) 



where A is a physical observable, — B is a perturbation 
potential, and is the Fourier-Laplace transform of the 
linear response function of A. (• ■ ■ ) F represents the av- 
erage in a NESS with a driving force (pump field) F and 
without the perturbation (probe force), and {»,»}p the 
Poisson bracket in classical systems (or the commutator 
divided by ih in quantum systems), which often reduces 
to an easily measurable quantity. 

In Ref. [12| they derived these properties assuming that 
a large system (which is composed of the system of inter- 
est, its environments and a driving source) is a (determin- 
istic) Hamiltonian system. The applicable range of the 
properties is quite large since almost no assumptions were 
imposed. It is not trivial, however, that the properties 
hold in stochastic systems. There are two points to ex- 
amine the validation in such systems. One is that explicit 
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derivation of the properties in stochastic systems extends 
the universality of the properties. The other is as follows. 
Some classes of stochastic systems may be regarded as ap- 
proximate models of certain large Hamiltonian systems, 
where the degrees of freedom in the environments are 
eliminated (although in some cases the connections be- 
tween these models are not sufficiently clear). Therefore 
the validation is utilized as a criterion for effectiveness of 
the approximate models because the properties should re- 
main valid if the approximations in reducing the original 
Hamiltonian models to the stochastic models are good. 

In this paper we show that the sum rule holds in 
the stochastic systems described by the Langevin equa- 
tions. It is valid in highly nonequilibrium steady states 
as well as in equilibrium states. Here, 'highly nonequilib- 
rium' is used in the following two senses: (1) the driving 
force may be arbitrarily large, and (2) the intensity of 
the Langevin noise is allowed not to satisfy the second 
fluctuation-dissipation relation (as in the case of nonequi- 
librium process in laser physics [H[ ). We also derive the 
asymptotic behavior, another general property shown in 
Hamiltonian systems [l2j], of response function in the 
Langevin systems. Furthermore we show that the sum 
rule for the velocity response function is violated in the 
overdamped cases. This implies that the overdamped 
model sometimes gives results inconsistent with those in 
larger Hamiltonian systems especially when they are re- 
lated to small-time-scale phenomena (which are correctly 
described by the Hamiltonian systems). We also discuss 
a relation between the sum rule and the Harada- Sasa 
equality 0, 0] , which is an expression of the energy dis- 
sipation in the Langevin models. 



II. MODEL 

We define a general form of the Langevin model. The 
examples are shown in Sec. IIVI 

We consider a system the state of which is speci- 
fied by a set of n stochastic variables denoted by C = 
(Ci) C2, • • • j Cn) (called system variables hereafter). The 
dynamics of the system is assumed to be characterized 
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by the following stochastic differential equation [l4| of a 
general form: 

n 

dQ(t) = M l (C(t);F)dt + Y,^AC(t);F) ■ dWj(t) 

i=i 

+ sff(t)K l (C(t))dt. (2) 

The first two terms in the right-hand side describe the un- 
perturbed dynamics of the system. F represents a pump 
field, which determines the degree of nonequilibrium of 
the system (e.g., an external driving force or a temper- 
ature difference). The first term is a deterministic part 
(Mi is a certain known function of £ and F). The second 
one represents the noise term, where Wi(t) is a Wiener 
process. We assume that the mean of Wi(t) is zero for 
any i and that Wi(t) and Wj(t) are uncorrelated if z ^ j. 
The symbol " • " implies the multiplication in the sense 
of Ito. The noise intensity is determined by Afij (which is 
a known function of £ and F) . The third term represents 
the probe force (perturbation) the response to which is 
of our interest (ff and i*Q are known functions of t and 
of £, respectively). We assume that in the unperturbed 
system (i.e., in the case of e = 0) is realized a certain 
steady state which is stable against perturbations. 

A situation to be supposed is as follows. The system 
of our interest is driven to a NESS by the pump field 
F. The steady state may be in a nonlinear response 
regime, in a linear response regime, or an equilibrium 
state, because F is allowed to be arbitrary (may be very 
large) and the noise intensity is arbitrary. To measure the 
response of the NESS, the probe force sffKj is applied 
to the system in addition to F. (This is sometimes called 
pump-probe experiment.) The linear response function 
to the perturbation effKj is defined by 



6(A(C(t))) 



eSff(s) 



(3) 



Here A is a physical observable which is a function of the 
system variables C, and y "/ Fe represents the average 
in a state with the parameter F and the perturbation. It 
should be noted that the response function thus defined 
is that of a NESS which may be in a nonlinear response 
regime because F is allowed to be large (irrespective of e) . 
[The meaning of "linear" in the linear response function 
of the NESS is the linear order in efjKj (not F).] 



III. GENERAL RESULT 

A. Sum rule for response function 

We now show the sum rule for the response function 
in the Langevin system. For this purpose we consider 
the stochastic differential equation for the observable A. 



This is given by the Ito formula [14j: 
<L4(C(*)) =AA(C(t))dt 



dA 



+ EE^i(C(*);^)^(C(*))-dw i (t) 



i=l 3 = 1 



OA 



•sEAW(c(*%(cw) d *. ( 4 ) 



1=1 



where the backward operator A is defined by 
n ft 

a = J2m<(C;F) ° 



d 2 



(5) 



i=lj=lfc=l 

Employing a method similar to that in Ref. @, we 
transform this differential equation into an integral form: 

A(C(t)) = e^ A A(C(t )) 

+ EE / udc(nF)^-^A(c(t>)) ■ dw 3 (t>) 



i=l j=l Jt ° 



*E/ fI(t')MC(t'))^re^ A A(C(t>))dt>, (6) 
l=1 J to °V 

where to is an initial time. For completeness we provide 
the detail of the transformation in the Appendix. 

The average of the second term in the right-hand 
side of Eq. (j6|) vanishes because the integrand is a non- 
anticipating function (the mean value formula of the Ito 
stochastic integral) [14] . Taking the average of Eq. ([6]), 
we thus obtain 

(A(C(t))) Fe = e^ A A(0 o )) 



-l J to 



+^E/ fmiMat'))^'^^)) dt ' 



da 



(7) 

By a functional differentiation of Eq. ([7]) with respect to 
ff (s), an expression of the response function is derived: 

< (t ~s;F) = (C(«)) ^e^ A A(C(s)) ' 

(8) 

for t > s (> t ), and (t - s; F) = otherwise. Note 
that the Langevin system satisfies the causality condi- 
tion. If s 3> to, this expression does not depend on s 
but on t — s only, because a steady state is realized in 
the unperturbed system. Thus we finally obtain the sum 
rule: 

Re*£>;i0— = $1(+0;F) 



*i(C)|£(C) 

°M / F,0 



(9) 
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where F) is the Fourier-Laplace trans- 

form of <£>^-.(t;F) and where wc have used 
f^Im&x^Lj; F)duj = (this is because $£.(r;F) is 
a real- valued function). This is the main result of the 
present paper. 

B. Asymptotic behavior 

In Ref. [Hj], in addition to the sum rule, an asymp- 
totic behavior of ^j^.(uj\F) was derived. We here show 
another derivation of it using the sum rule. We first note 
that the response function satisfies a dispersion relation 
if the system satisfies the causality condition; 

/OO -1 J / 
R^Uu';F)— , 
-OO * 

(10) 

where denotes the principal value. By multiplying oj to 
the both sides of this equation and by taking the u> — > oo 
limit, we obtain the asymptotic behavior: 

/OO J / 

Rel>£ (u/;F) — 
-oo ' 7T 

where we have exchanged the integration and the limit 
procedure in the first equality and have used the sum rule 
in the second equality. Hence the asymptotic behavior 
is valid in the systems (including the Langevin systems) 
where the sum rule and the dispersion relation hold. The 
asymptotic value is the same as the sum value in the sum 
rule. 

C. Remarks 

Here we make remarks on the results. 

First, Eq. © should be interpreted as a prediction 
on the sum value (integral value) of $^-.(w;F) at many 
different values of uj [therefore Eq. © is called sum 
rule]. $k (w;F) for each w is easily measured by ex- 
periments or numerical simulations, e.g., in the following 
way: apply a sinusoidal probe force with frequency oj [i.e., 
f p (t) = sinwt], measure the observable A for sufficiently 
long time (^ 1/w), and calculate the Fourier component 
(Aui) F,e at oj from the time-series data of A ($^. (oj; F) — 

lim e ^o(^4-cj)_F,e/£)- The sum rule [Eq. (j9))] states that the 
sum value measured in the system with a small perturba- 
tion equals the average value of 5(C) = Kj(C,)dA(C,)/dC,j 
measured in the system without the perturbation. In 
many cases 5 is a quantity that is easy to measure. 
Therefore it is possible to test the sum rule in experi- 
ments and simulations. (In Ref. [l2| we demonstrated 
the validity of a sum rule by a numerical simulation.) 



Equation (|9]) also seems to be a prediction on the value 
of the response function immediately after the probe 
force is applied, as seen in the middle of Eq. In 
contrast to the measurement of 5>^. (ui; F), however, it is 
very difficult to measure 3?^, (+0; F) in experiments and 
numerical simulations. Therefore it is more natural to 
interpret Eq. (JSJ as a sum rule for $^ (oj; F). 

Second, the sum rule holds for any steady states (if 
they are stable) including the states in a nonlinear re- 
sponse regime (at large F) as well as an equilibrium state 
(F = 0). Although the form of the rule is the same for 
all the states, the value of the sum may vary as F is 
changed. This is because the steady state distribution 
function, which appears in averaging 5(C) m Eq. ©, is 
dependent on F. 

Third, the expression of the response function [Eq. ([8])] 
in more concrete examples was derived in some literatures 
HE]]) although they did not derive the sum rule from it. 
The expression itself is not so useful, because it is hard 
to calculate in analytic way unless one knows the ex- 
plicit form of the steady state distribution function, and 
because it is hard to measure in experiments and numer- 
ical simulations due to the complicated factor e^ s ^ A . In 
contrast, the rightmost side of Eq. © is much easier to 
measure since it does not contain such factors, and tests 
of the validity of Eq. © are possible as mentioned in the 
first remark. 

Fourth, the sum rule is different from the moment 
sum rule (e.g., Ref. Q). The statement of the mo- 
ment sum rule is that the sum value of (oj;F) [and 
lu x <&k. (oj; F)} converges to a certain value. However, the 
dependence of this value on F and the other parameters 
is not given by it, although the dependence on F is most 
interesting point in nonequilibrium statistical mechanics. 
In contrast, the statement of the sum rule is not only the 
convergence of the sum value but also its equivalence to 
the average value of 5(C), which gives the dependence of 
the sum on several parameters (including F). 

5 is independent of the state and the system (without 
probe). (It depends on the probe force and the observ- 
able of our interest, both of which may be chosen irre- 
spective of the state and the system.) Therefore there 
are no differences in 5 between equilibrium states and 
NESSs and between non-interacting systems and inter- 
acting systems. The dependence on the state and the 
system appears only in averaging it, which leads to the 
dependence of the sum on F in general. When 5 is inde- 
pendent of C: in particular, this dependence disappears 
and the sums are the same for any steady states and 
for any systems. This fact cannot be predicted by the 
moment sum rule. 



IV. EXAMPLES 

In this section we describe some examples of the 
Langevin model defined in the general form [Eq. ^] in 
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Sec. HH and see the concrete forms of the sum rule in the 
examples. It should be noted that in Eq. (J2J the relevant 
quantities to the sum rule are K^s. 

Furthermore we show that the sum rule for the velocity 
response function is violated in the overdampcd case. 



A. Underdamped case 

One of the simplest examples is the single-particle un- 
derdamped Langevin model in one dimension which is 
described by the following equations: 



d P (t) 
dt 



-v{t) 

m 



dx(t) p(t) 



dt 

Here, p, x, 



(12) 
(13) 



m, 



and 7 are the momentum, position, 
mass and friction coefficient of the particle, respectively. 
U(x) is a potential, F is an external driving force, and 
—B(x) is a perturbation potential of the probe. £(t) 
is a white Gaussian noise with zero mean and satisfies 
(£(*)£(*')) = 2DS(t - t'). F is allowed to be arbitrarily 
strong. A nonequilibrium steady state in the nonlinear 
response regime is realized for large F, while an equilib- 
rium state is realized for F = 0. 

More precise forms of Eqs. (fT2")l and (|13l) are given by 
Eq. ([2J where n = 2, £i = p, and £2 = x, and where 

M 1 (p,x;F) = -±p- — (x) + F, M 2 (p, x; F) 



Afu(p,x;F) 



m 



dx 

K 2 (p,x) = 0. 



P_ 

m ' 



0, 



In this case the sum rule reads 



Re$^ 



dio 



8B dA s 



(14) 



F,0 



This form is the same as that in the classical Hamiltonian 
models where the systems of interest are single-particle 
systems. In particular, when considering the momentum 
response to the spatially homogeneous probe force (i.e., 
A = p and B — x), we have the F -independent sum 
value: 

Re*£(w;F)— = 1. (15) 

The second example is a three-dimensional many- 
particle underdamped Langevin model described by 

d&(t) 



dt 



dU ,_ , .v dV . . ,v 

^(^W)-^r(KW}) 



dt 



+ F + Ut) 



f) R 

eF{t) — (p,{tlr,(t)), 



.(*)), (16) 
(17) 



for = 1,2, ••• , iV. Here, a = (a 1 , a 2 , a 3 ) represents a 
three-dimensional vector, p^, r^, m M , and 7^ are the 
momentum, position, mass and friction coefficient of the 
/j,th particle, respectively. U(r) is a single-particle po- 
tential and V({7v}) is an inter-particle potential. F is 
an external driving force, the strength of which is arbi- 
trary. —B{p,r) is a probe potential (perturbation) which 
is assumed to depend on a momentum as well as on a 
position. (An example is an interaction of an electron 
with an electromagnetic field.) is a white Gaus- 

sian noise with zero mean and satisfies (£')) = 

2D5 aa 'S t2t2 >S(t — t') (a, a' = 1,2,3). More precise forms 
of these equations are given by Eq. © with n = 6N, 
Ce^-e+a = p%, and = r% {a = 1,2,3). The 

relevant quantities AVs to the sum rule are written as 

K 6l/t - 3+a ({p u ,f v }) = J^(p M ,^). 

We thus obtain the sum rule for the response function 
$g(t- s;F) = S(A{{p,(t),r,(t)})) FQ /e5p( s )\ £=0 



'hb^( W ;F) — 

7T 



dr 

\ fJ,,Ot u " 



(18) 



F,0 



This form is the same as that in the classical Hamiltonian 
models [Eq. flTJ]. This reduces to Eq. (jUJ) if N = 1, 
dB I dp = and the motion is restricted to one dimension. 

In Langevin models the noise intensity D is usually as- 
sumed to be equal to ^k^T [which is called the second 
fluctuation-dissipation relation (FDR) [l|], where T is the 
temperature of the environment and k-Q is the Boltzmann 
constant. However, as seen in the general result and in 
the above two examples, this assumption is not neces- 
sary for the validity of the sum rule because the noise 
intensity does not explicitly contribute to the rule. This 
indicates that the sum rule holds in a very wide range of 
systems from highly nonequilibrium systems (in the sense 
that the second FDR is violated) to purely deterministic 
systems (D = 0). The violation of the second FDR is of- 
ten seen in the treatment of nonequilibrium process (far 
from equilibrium) in light-emitting devices by quantum 
Langevin equation (l3| . 

For the same reason the noise intensity D may depend 
on F. Therefore the treatment in the present paper is 
applicable also to the systems in which the second FDR 
is gradually violated as F increases. Also D may depend 
on the position of the particles in the two examples and 
on /1 and a in the second example. One of the conse- 
quences from this fact is that the sum rule is valid also in 
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heat conducting nonequilibrium systems driven by tem- 
perature difference. It should be noted that the valid- 
ity of the sum rule is independent of the noise intensity 
whereas the value of the sum depends on it because the 
steady states depend on it in general. 

Finally we make a brief comment. In some Langevin 
systems the natural interpretation of the Langevin equa- 
tions used in physics is the Stratonovich-type stochas- 
tic differential equations [ll|. In such cases we should 
first interpret the given Langevin equations to the cor- 
responding Stratonovich-type equations, then transform 
them into the Ito-type equations [Eq. ||2J|], and finally ap- 
ply the general form of the sum rule. Note that the form 
of the sum rule is identical irrespective of the interpre- 
tations, although the value of the sum depends on them 
because the steady states may be different in different 
interpretations. In the above examples (with position- 
independent D), the values as well as the forms of the 
sum rules are independent of the interpretations because 
the forms of the Stratonovich-type and Ito-type equa- 
tions are the same except for the senses of the multipli- 
cations in these cases. 



B. Overdamped case 

In this subsection we discuss the sum rule in the over- 
damped case. For simplicity we consider a single-particle 
model (the extension to many-particle models is straight- 
forward). We consider the single-particle overdamped 
Langevin model in one dimension which is described by 



7" 



dx(t) 
~dT~ 



dU 
dx 



(x(t))+F + Z(t) 



(19) 



Here the notations are the same as those in Eqs. ([12)) and 
(|13[) . £(t) is a white Gaussian noise with zero mean and 
satisfies (£{t)£(t')) = 2DS(t — t'). A more precise form of 
Eq. (Tl9l) is given by Eq. ([2} where n = 1 and Ci = x, and 
where 



Mi (x; F) = -- 

7 



7 

1 dB 
K x (x) = -j-ix). 
7 ox 

Even in this model the sum rule is valid if the observable 
A is a function of only x; 



du 



1 /dB, s dA 



Re^F)- = -( K - { x)-( X ) 



(20) 



F,0 



One sometimes considers as an observable of interest 
the velocity v, which is defined by v(t)dt — dx(t) in the 
overdamped model. Because we have not derived the 
sum rule in the cases that the observable is a function of 



the differentials of the system variables, the general result 
[Eq. ©] does not ensure that the sum rule holds for the 
response function & B of the velocity in this model. How- 
ever, from the viewpoint that the overdamped model is 
regarded as a coarse-grained model of the underdamped 
one, <fr v B must satisfy the sum rule if the coarse-graining 
procedure is good. 

To investigate this point we again consider the under- 
damped model described by Eqs. (fT2[) and (|T3|) . Because 
the momentum of the particle is included in the system 
variables in the underdamped model, we can safely con- 
sider the sum rule for the response function of the velocity 
p/m, which reads 



Re$ 



p/j 



_.dw 7 1 /dB 

(ui;F) — = ( — 

7r mj\dx 



(21) 



F,0 



The right-hand side diverges in the overdamped limit 
(m/7 — > 0). Therefore a necessary condition for the va- 
lidity of the sum rule for $^ is that $^(+0; F) is diver- 
gent in the overdamped model. We next examine this 
condition directly in the overdamped model. § V B is cal- 
culated by averaging Eq. flT9l and by functionally differ- 
entiating the result with respect to f(s): 



5 I dx(t) \ 



S]F) eSf(s){ di'/j 



1 /dB 



7^ \ dx 
1 /dB 



F.O 



(22) 



where we have used the expression [Eq. ©] in the first 
term in the rightmost side. Owing to the second term, 
$ V B (+0;F) is divergent [l6j |. Thus the sum rule for <& V B 
seems to be valid in the overdamped model in the sense 
that the both sides of it are divergent. More careful anal- 
ysis, however, reveals that the diverging behaviors are 
different between Eq. ([2"2"[) and the overdamped limit of 
Eq. (f2T|) . The diverging behavior of Eq. ([22} is domi- 
nated by the delta function in the second term in the 
rightmost side. This comes from the functional derivative 
Sf(t)/8f(s), the order of which is estimated as 0(1/ At). 
Here At is the smallest time scale of our observation on 
the system. On the other hand, the diverging behav- 
ior of the overdamped limit of Eq. ([211) is dominated by 
the factor 7/771 in the right-hand side. Since the over- 
damped Langevin model should be interpreted as an ef- 
fective description of the underdamped one in the time 
scale At ^> m/7, one must first take the m/j — > limit 
and then take the At — > limit to have continuous time 
limit of the overdamped model. Therefore the diverging 
behavior is stronger in the overdamped limit of Eq. ([21[) 
[0(7/ra)] than in Eq. ||22J) [0(1/ At)]. In this sense the 
sum rule for & B is violated in the overdamped model. 
This is consistent with the fact that the overdamped 
Langevin model is valid (as a coarse-grained model of 
the underdamped one) only in the frequency range of 
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w < 7/m (which would result in an incorrect contribu- 
tion to the sum from the higher frequency region). 



V. RELATION TO THE HARADA-SASA 
EQUALITY 

In Refs. @, Harada and Sasa derived an equality 
(the Harada-Sasa equality) on the energy dissipation rate 
in nonequilibrium Langevin systems. We here discuss the 
relation between this equality and the sum rule. 

For simplicity we consider the single-particle under- 
damped Langevin model described by Eqs. (fT2|) and (fT3| . 
In this model the energy dissipation rate J from the sys- 
tem to the environment is defined by [l5|, HH 



j(t)dt = (i P (t) - e(<)) o dx(t), 



(23) 



where the symbol " o " represents the multiplication in 
the sense of Stratonovich. The Harada-Sasa equality is 
an expression of the average of J: 



7 

(J)f,o = — 9 (p)f,o 



C p (w; F) - 2mk B TRe^(u}; F) 



du> 

2tt' 

(24) 



where the usual assumption (the second FDR) D = ^k^T 
on the noise intensity is imposed. C p (uj; F) is the Fourier 
transform of the time-correlation function C p (t;F) of 
the momentum; C p {t;F) = ((p(r) - (p) F ,o)(p(0) ~ 

(P>F,o)> F)0 - 

We show that this equality is derived with the help 
of the sum rule. First we note that Eq. (|23|) is rewrit- 
ten in the Ito type: J(t)dt = (jp 2 {t)/m 2 — D/m)dt — 
(y/2D/m)p(t) ■ dW(t). This is derived by substituting 
"£(*) = \/2DdW(t)/dt" and dx(t) = p(t)dt/m into 
Eq. (|23l) . and by using the Stratonovich- Ito transforma- 
tion [14{ . Since the average of the last term vanishes due 
to the mean value formula, the average dissipation rate 
is written in a simple form [l5|, [l9| : 



(J)f,o 



;(P 2 



F,0 



D 

rn 



(25) 



By multiplying Re$P(w; F)du}/n, which is equal to 
1 owing to the sum rule [see Eq. (fT5|) ] . to the last term, 
and by noting (p 2 ) F ,o = (p) Ffl 2 + CP (w; F)6u/2ir, 
we obtain 



+ 



J_ 

rn 2 

7 

m 2 



(p)f/ 



2mD 

7 



Re<b p (u;F) 



dco 
27r" 
(26) 



This becomes the Harada-Sasa equality if D = ^U-qT [2C 



As seen in the above derivation, the validation of the 
Harada-Sasa equality requires that Eqs. (H"5l and (|2"5)l 
hold. The former holds in a wide range of nonequilibrium 
systems (even in systems other than Langevin systems) 
because it is a specific form of the sum rule. In con- 
trast, the validity range of the latter is not so large. The 
meaning of Eq. (|25p with D = jk-gT is that the ratio of 
the average dissipation rate and the difference between 
the kinetic energy of the system and the environment 
temperature, {J) f,o/ ({p 2 ) Fft/m — AtbTJ, is a constant, 
7/m, which is independent of F. In some systems other 
than Langevin systems, however, this does not hold, es- 
pecially in states far from equilibrium. For example, in a 
numerical simulation of a model of electrical conduction 
[Ml2l, it is clearlyseen that (J) F ,o/((p 2 )F,o/m- k B T) 
does depend on F [231 ] . In this sense the validity range of 
the Harada-Sasa equality is restricted to the systems in 
which the ratio is constant. A sufficient condition for this 
might be the distinct separation of time scales, as they 
mentioned in Ref. Q. It should be noted, however, that 
even when the ratio is not constant there remains an- 
other possibility. That is, it might be possible that one 
can define an independent j(F) in a certain way and 
that "f{F)/m equals the ratio (J)f,o/ ((p 2 )-F,o/m — 
in the steady state at each F. Whether this is true or 
not should be tested in systems (e.g., in the model in 
Refs. [U Ull) which cannot be described by Langevin 
equations. 



VI. SUMMARY 

In this paper we extended the validity range of the sum 
rule (and the asymptotic behavior) for the linear response 
function of steady states to a class of stochastic models 
described by a general form of the Langevin equations. 
This holds for a wide range of the steady states (if they 
are stable) including highly nonequilibrium states as well 
as equilibrium states, because the driving force is allowed 
to be large (e.g., to be in a nonlinear response regime) 
and the noise intensity may be arbitrary (e.g., not to sat- 
isfy the second fluctuation-dissipation relation) . The sum 
rule is a property which normal nonequilibrium models 
should have and therefore is used as a touchstone to ex- 
amine the correctness of results in experiments and the- 
ories of NESSs. In the overdamped Langevin model the 
sum rule for the velocity response function does not hold, 
which suggests that results in the overdamped model 
should be treated with care if they are concerned with 
small time scales. 

We also showed the relation of the sum rule to the 
Harada-Sasa equality. The equality is reduced to a sim- 
pler form when one uses the sum rule with a specific 
choice of observable. 

Further extension of the validity range of the sum 
rule remains as theoretical issues. Extension to time- 
dependent case where Mj and Mij in Eq. explicitly 
depend on t is straightforward. It would be also in- 
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teresting to examine the validation of the sum rule in 
non-Markovian models. It is easily extended to the non- 
Markovian models which become Markovian if appropri- 
ate dynamical variables are added to the original system 
variables. For other non-Markovian models, the method 
used in Refs. [3, [HI to generalize the Harada-Sasa equal- 
ity to non-Markovian cases might give hints. Finally, it 
is also important to investigate higher-order responses of 
NESSs in stochastic systems as studied in Hamiltonian 
systems recently [26j . 
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Appendix: Derivation of Eq. ([6]) 



Because A is arbitrary the above equation is regarded as 
a stochastic differential equation for fl. In order to trans- 
form this equation into an integral equation, we introduce 
an operator, Q[t) = n(t\e~^~ to ^. From Eq. (1X11) the 
differential equation for Q is given by 



i=i j=i 



-(t-to)A 



-(t-to)A 



dt. (A.2) 



Then formally integrating this equation from to to t with 
the initial condition H(to) = 1 and multiplying e^~ to ^ A 
from the right, we obtain 



We introduce the time-evolution operator Q(t) corre- 
sponding to Eq. for an arbitrary function C of (, the 
value at C(t) is given by C(((t)) = f2(t)C(((t )), where 
to is the initial time and i2(to) = 1. Then Eq. §4$ is 
rewritten as 



df2(t)A(C(t )) = m [AA\ (C(t ))dt 

dA 



i=i j=i 



dA 



{C{to))-dW 3 (t) 



(C(h))dt. 



Q{t) = e {t - to)A 

+EE / 0(0 



i=i j=i ut « 

eE / f!(t'W> 



dt'. (A.3) 



(A.l) Acting this equation on A(((to)), we finally get Eq. ©. 
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